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I. INTRODUCTION 



o 
o 
o 

Various mechanical problems can be elegantly approached by the Hamiltonian formalism, which not only found 
• well-established ground in classical theories , but also provided much physical insight in the early development of 
[ quantum theories It is curious though that the concept of canonical transformations, which plays a fundamental 

■ role in the Hamiltonian formulation of classical mechanics, has not attracted as much attention in the corresponding 
I formulation of quantum mechanics. A relatively small quantity of literature is available as of now on this subject 

[4-10]. The main reason for this is probably that canonical variables in quantum mechanics are not c- numbers but 
CNJ noncommuting operators, manipulation of which is considerably involved. In spite of this difficulty, the great success 
^ of canonical transformations in classical mechanics makes it desirable to investigate the possibility of application of 
the concept of canonical transformations in quantum mechanics at least to the extent allowed in view of the analogy 
with the classical case. 

The usefulness of the classical canonical transformations is most visible in the Hamilton-Jacobi theory where one 
] seeks a generating function that makes the transformed Hamiltonian become identically zero [Q . A quantum analog 
^\ of the Hamilton-Jacobi theory has previously been considered by Leacock and Padgett |^ with particular emphasis 
0^ ' on the quantum Hamilton's characteristic function and applied to the definition of the quantum action variable and 
, the determination of the bound-state energy levels ||ll[]. However, the dynamical aspect of the quantum Hamilton- 
^ ' Jacobi theory appears to remain untouched. In the present study, we concentrate on this aspect of the problem, 
^ I , , and derive the time-dependent quantum Hamilton-Jacobi equation following closely the procedure that lead to the 
(~| ' classical Hamilton-Jacobi equation. 

[ The analogy between the classical and quantum Hamilton-Jacobi theories can be best exploited by employing 
^ ^- the idea of the quantum generating function that was first introduced by Jordan Q and Dirac and recently 
reconsidered by Lee and I'Yi The "well-ordered" operator counterpart of the quantum generating function is 

used in constructing our quantum Hamilton-Jacobi equation, which resembles in form the classical Hamilton-Jacobi 
equation. By means of well-ordering, a unique operator is associated with a given c-number function, thereby the 
ambiguity in the ordering problem is removed. We identify the quantum generating function accompanying the 

■ quantum Hamilton-Jacobi theory as the quantum Hamilton's principal function, and apply this theory to find the 
dynamical solutions of quantum problems. 

The prevailing conventional belief that physical observables should be Hermitian operators invokes in our discussion 
the unitary transformation that transforms one Hermitian operator to another. This along with the fact that the 
unitary transformation preserves the fundamental quantum condition for the new canonical variables [Q, P] — ih 
if the old canonical variables satisfy [q,p] = ih provides a good reason why we call the unitary transformation the 
quantum canonical transformation. This definition of the quantum canonical transformation is analogous to the 
classical statement that the classical canonical transformation keeps the Poisson brackets invariant, i.e., [Q,P]pb — 
[q,p\pB = 1. In our current discussion of the quantum canonical transformation we will consider exclusively the case 
of the unitary transformation. 

The paper is organized as follows. In Sec. II the quantum canonical transformation using the idea of the quantum 
generating function is briefly reviewed, and the transformation relation between the new Hamiltonian and the old 
Hamiltonian expressed in terms of the quantum generating function is derived. From this relation, and by analogy 
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with the classical case, we arrive at the quantum Hamilton- Jacobi equation in Sec. III. It will be found that the 
unitary transformation of the special type U{t) — T{t)A where T{t) is the time-evolution operator and A is an 
arbitrary time-independent unitary operator satisfies the quantum Hamilton-Jacobi equation. Sec. IV is devoted to 
the discussion of the quantum phase-space distribution function under canonical transformations. The differences 
between our approach and that of Ref. [ p^ are described. Boundary conditions and simple applications of the theory 
are given in Sec. V, where to perceive the main idea easily most of the discussion is developed with the simple case 
A = I, the unit operator, while keeping in mind that the present formalism is not restricted to this case. Finally, Sec. 
VI presents concluding remarks. 



II. QUANTUM CANONICAL TRANSFORMATIONS 

Let us begin our discussion by reviewing the theory of the quantum canonical transformations ||5|,[l0|]. A quantum 
generating function that is analogous to a classical generating function is defined in terms of the matrix elements of 
a unitary operator as follows j|] , 

e^^^^'^^^Q^^'^'^ ^ {qi\Q2)t ^ {qi\U{t)\q2), (1) 

where the unitary operator U{t) transforms an eigenvector of q into an eigenvector of Q = UqU\ i.e., \Qi)t = U{t)\qi) 
(and |i-i)f = ?7(t)|pi)).|^ Different types of the quantum generating function can be defined similarly i.e., 

e«F.to,p.,t)A = (g^|P2), ^ {q^\U{t)\p2), e'^^^P^'Q-'^y^' = (pi|Q2>t = (pi|;7(i)|<72>, and eP^^vM)!^^ ^ (p^jp^), ^ 
(Pilf^Wb). 

The quantum canonical transformation, or the unitary transformation, corresponds to a change of representation or 
equivalently to a rotation of axes in the Hilbert space. The unitary transformation guarantees that the fundamental 
quantum condition = [9,p] = ih holds, the new canonical variables {Q,P) are Hermitian operators, and the 

eigenvectors of Q or P form a complete basis. One should keep in mind that the eigenvalue Qi has the same numerical 
value as the eigenvalue qi because the unitary transformation preserves the eigenvalue spectrum of an operator 
In cases where it is convenient, one is free to interchange qi (pi) with Qi (Pi). 

Transformation relations between {q,p) and {Q,P) can be expressed in terms of the "well-ordered" generating 
operator Fi{q, Q, t) |^ that is an operator counterpart of the quantum generating function Pi(<?i, Q2, i) as follows:|^ 

. 9Pi(g,Q,t) ~ dFi{q,Q,t) 

dq dQ 

Similar expressions for other types of the generating operators can be immediately inferred by analogy with the 
classical relations. For a later reference, we present the relations for P2('?, P, t) below, 

. dF2{q,P,t) ~ dF2{q,P,t) 

dq dP 

It is interesting to note that, whereas the four types of the generating functions in classical mechanics are related 
with each other through the Legendre transformations |^] , the relations between the quantum generating functions of 
different types can be expressed by means of the Fourier transformations. For example, the transition from Fi{qi,Q2, t) 
to P2(<Zi, P2, t) can be accomplished by 

e'P-<^''^^P-''y^ ^ J dQ2{qi\Q2)tt{Q2\P2)t, 



^ An eigenvalue Xi and an eigenvector \Xi) of an operator X are defined by the equation, X\Xi) — Xi\Xi) {X = q, p, Q, 
and P). Different subindices are used to distinguish different eigenvalues or eigenvectors, e.g., X2, \X2), etc.; The subscript t 
on a ket |)t (bra t{\) expresses time dependence of the ket |)t (bra t{j). 

^ A well-ordered operator G{X,Y) is developed from a c-number function G{Xi,Y2) such that {Xi\G{X ,Y)\Y2) = 
GiXi,Y2){Xi\Y2} §. For example, ff G(Xi, Y2) = X1Y2 + Y^Xf, then G{X, Y) = XY + X^Y^. 
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The usefulness of the concept of the quantum generating function can be revealed, for example, by considering 
the unitary transformation U = e^^^^^^^ where g is an arbitrary real function. From the definition of the quantum 
generating function, we have 



The well-ordered generating operator is then given by 



F2iq,P)=g{q) + qP + i^ln2wh, (6) 
and Eq. (3) yields the transformation relations 

Q = q, (7) 

P = (8) 

This shows that, in some cases, an introduction of the quantum^ generating function can provide an effective method 
of finding the transformation relations between and {Q,P) without recourse to the equations Q = UqU^ and 
P = UpU''. 

Now we consider the dynamical equations governing the time-evolution of quantum systems. The time-dependent 
Schrodinger equation for the system with the Hamiltonian H{q,p, t) is given in terms of a time-dependent ket \tp)t by 

in^^\lb)t=H{q,p,tmt. (9) 
In Q-representation the time-dependent Schrodinger equation takes the form 

th-^^i'Q{Qi,t) = K (Q,,-rn-^^,?j V'3(Qi,t), (10) 
where tp'^{Qi,t) = t{Qi\tp)t, and 

K{Q,P,t) = H{q,p,t)+mU^. (11) 

The second term on the right hand side of Eq. (11) arises from the fact that we allow the time dependence of the 
unitary operator U{t), which indicates that, even though we adopt here the Schrodinger picture where the time 

dependence associated with the dynamical evolution of a system is attributed solely to the ket |V')t, Q and \Qi)t may 
depend on time also. In terms of the generating operator Fi{q, Q, t), Eq. (11) can be written as 

K{Q, P, t) = H{q, p, t) + ^IlMA . (12) 

The equivalence of Eqs. (11) and (12) can be proved as shown in Appendix A. It is important to note that K{Q, P, t) 
plays the role of the transformed Hamiltonian governing the time-evolution of the system in Q-representation. The 
analogy with the classical theory is remarkable. 

III. QUANTUM HAMILTON-JACOBI THEORY 

We are now ready to proceed to formulate the quantum Hamilton-Jacobi theory. One can immediately notice 

that, if K{Q,P,t) of Eq. (12) vanishes, the time-dependent Schrodinger equation in Q-representation yields a simple 
solution, ijj'^ = const. This observation along with Eq. (2) naturally leads us to the following quantum Hamilton- 
Jabobi equation. 
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where, following the classical notational convention, we denote the generating operator that is analogous to the classical 
Hamilton's principal function by Si{q, Q, t). Eq. (13) bears a close formal resemblance to the classical Hamilton- Jacobi 
equation. It, however, differs from the classical equation in that it is an operator partial differential equation. The 
procedure of solving dynamical problems is completed if we express the wave function in the original g-representation 
as 

= J {qi\Q2)tt{Q2mtdQ2, 

= y e'^^ V^(Q2)rf(32, (14) 

where Si{qi, Q2,t) is the c- number counterpart of Si{q,Q,t), and is obtained by replacing the well-ordered q and 
Q in Si, respectively, with qi and Q2- In Eq. (14), t is dropped from ip'^ , since t{Q2\'ip)t = const. As is the case 
for the classical Hamilton-Jacobi equation, the mission of solving dynamical problems is assigned to the quantum 
Hamilton- Jacobi equation. 

Even though we arrive at the correct form of the quantum Hamilton-Jacobi equation, it seems at first sight quite 
difficult to attain solutions of it due to its unfamiliar appearance as an operator partial differential equation. Thus 
it seems desirable to search a corresponding c- number form of the quantum Hainilton- Jacobi equation.^ For this 
task, we note that, if the unitary operator U{t) is assumed to be separable into U{t) = T{t)A, where T{t) is the 
time-evolution operator and A is an arbitrary time-independent unitary operator, then il)^{Qi,t) = t{Qi\'4')t = 
{qi\A'^f^{t)f{t)\tJj{t = 0)) = ((?i|it|Vj(f = 0)) = const. This means that the left hand side of Eq. (10) becomes zero, 
i.e., the canonical transformation mediated by a separable unitary operator is exactly the one that we seek. Assuming 
U{t) = f{t)A, we rewrite Eq. (1) as 

giS^iq^,Q2,t)/h ^ (^q^\f^t)A\q2). (15) 

Differentiating this equation with respect to time, wo obtain 

Eq. (16) leads immediately to the desired c-number form of the quantum Hamilton-Jacobi equation 



fTl -t, ^ dSi{qi,Q2,t) 
H qi,-in- — ,t + 



dqi J dt 



gjSi(gi,Q2,t)/ft _ Q_ (^ij-^ 



Substitution of S2{qi, P2,t) for 5*1(51, Q2, t) generates another c-number form of the quantum Hamilton-Jacobi equa- 
tion. The equations for the cases of ^3(^1, Q2, t) and Si{p\,P2,t) can be derived through a similar process. 
Consider a one-dimensional nonrelativistic quantum system whose Hamiltonian is given by 

H{q,p,t) = ?^ + V{q,t). (18) 
The c-number form of the quantum Hamilton-Jacobi equation (17) for this problem becomes 



o(^) -^o^ + nai,i) + ^ = o. (19) 



We can see clearly that, in the limit h ^ 0, the above equation reduces to the classical Hamilton-Jacobi equation. 
The second term of Eq. (19) represents the quantum effect. We note that it has been known from the early days that 
substitution of ip{q,t) — e^s{q,t)/h -j^^^ ^^le Schrodinger equation gives rise to the same Hamilton-Jacobi equation for 
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S(a, t)^ where S(q,t) is interpreted merely as the complex- valued phase of the wave function (see, for example, Ref. 
|13| ) . The present approach more clearly shows the strong analogy between the classical and quantum Hamilton- Jacobi 
theories emphasizing that the quantum Hamilton's principal function S'l which is related with the wave function via 
Eq. (14) plays the role of the quantum counterpart of the classical generating function. Moreover, as discussed later 
in Sec. V, e^^^l^ defined in Eq. (14) can be interpreted as a propagator under a certain choice of A. 

It may be viewed that the Hamilton- Jacobi equation in the form of Eq. (19) is no more tractable analytically than 
the Schrodinger equation for general potential problems. Nevertheless, it would be possible at least to obtain an 
approximate solution of it using a perturbative method as follows. Since the solution of Eq. (19) is given by the 
classical Hamilton's principal function in the limit ?i — > 0, we can expand the general solution in powers of h: 

5i = 5j"^-f?i5^'^+n2^f^+..., (20) 

where Sf^ is the classical Hamilton's principal function. Substituting Eq. (20) into Eq. (19) and collecting coefficients 
of the same orders in ?i, we can obtain 

1 (dSf^V , , dSf^ , , 

and 



2^ dqx dqx 2 dq\ dt 



= 0, n > 1. (22) 



Given the solution s["^ of the classical Hamilton-Jacobi equation (21), we solve Eq. (22) to find s[^\ s[^^ can be 
determined subsequently from the knowledge of s[^^ and s[^\ and so forth. We note that Eq. (22) is linear in S'j;"'' and 
first-order differential in qi for s["'\ Thus, from a practical viewpoint, Eqs. (21) and (22) could be more advantageous 
to deal with than Eq. (19) as long as the classical Hamilton's principal function that is the solution of Eq. (21) is 
readily available. 

The present formalism provides an encouraging point that the well-ordered operator counterpart of the quantum 
Hamilton's principal function gives also the solutions of the Heisenberg equations through Eq. (2). If we consider the 
case tj{t) — T{t), we can obtain in the Heisenberg picture the relations {(IHtPh) = {T^qsT,T^psT) and {Qh,Ph) = 
{f''Qsf,f''Psf) = {f'^fqsf'^f,f^fpsf'^f) = {qs,Ps), where we attached the subscript s and h to operators to 
explicitly denote, respectively, the Schrodinger and the Heisenberg pictures. Thus, when expressed in the Heisenberg 
picture Eq. (2) turns into 

dSi{qH,qs,t) . dSi{qH,qs,t) 

PH = , PS = , (23) 

oqH oqs 

and from these transformation relations we can obtain qn and pn as functions of time and the initial operators qs 
and pg. Obviously, qHiQs^PSj t) and pHi^s^PSi t) obtained in this way evolve according to the Heisenberg equations. 



IV. QUANTUM PHASE-SPACE DISTRIBUTION FUNCTIONS AND CANONICAL 

TRANSFORMATIONS 

Since our theory of the quantum canonical transformations is formulated with the canonical position q and momen- 
tum p variables on an equal footing, it would be relevant to consider the phase-space picture of quantum mechanics, 
exploiting the distribution functions in relation to the present theory. 



^For a stationary state of a system whose Hamiltonian does not depend explicitly on time, one may put S{q,t) = W{q) — Et 
and obtain a differential equation for W{q). To find a solution to the resulting equation, one may then use the expansion of W 
in powers of h. This approach has been extensively considered in connection with the well-known WKB approximation. In the 
present paper, the formalism is developed for general nonstationary states (of systems that can possibly have time-dependent 
Hamiltonians) . 
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A. Distribution functions 



For a given density operator p, a general way of defining quantum distribution functions proposed by Cohen is 
that 

F^qi^Put) J J dxdydg2(g2 +2/|p|g2 -2/)/(x,22//?i)e-fe-'^)e-*2^^'^/^ (24) 

Various choices of f{x,2y/h) lead to a wide class of quantum distribution functions [Tsl] . To mention only a 
few, the choice f — 1 produces the well-known Wigner distribution function ||l6|, while the choice f{x,2y/h) — 
g-hx /4ma-may jti yj^Ug ^^j^q Husimi distribution hmction that recently has found its application in nonlinear dynam- 
ical problems [0. The transformed distribution function is defined in (Qi, Pi) phase space likewise by 

G-^(gi,Pi,t) = ^y y J dXdYdQ2t{Q2+Y\p\Q,-Y)tf{X,2Y/h)e^''^'^--'^^h-''''''^/\ (25) 

Our main objective here is to find a relation between the old and the transformed distribution functions. After a 
straightforward algebra, which is displayed in Appendix B, it turns out that the transformation relation between the 
two distribution functions can be expressed as 

Gf{Ql,Pl,t)^ J J dq2dp2K{Qi,Pi,q2,P2,t)Ff{q2,P2,t), (26) 

where the kernel k is given by 

K{Qi,Pi,q2,P2,t)^-\- f f f f f f dXdYdQ2dxdydJ^^''^^^^^ 



^TT^nJ J J J J J " f{x,2y/h) 

This expression for the kernel can be further simplified if integrations in Eq. (27) can be performed with a specific 
choice of the function /. For instance, the simple choice f — 1 provides the following kernel for the Wigner distribution 
function, 

K{Qi,Pi,q2,p2,t)^^ [ /'drdye^[^i(?^-«'«i-^'*)-^i*te+^^'3i+^'*)le*[^P^-^^il. (28) 



This equation was first derived by Garcia-Calderon and Moshinsky |18| without employing the idea of the quantum 
generating function. Curtright et al. p9| also obtained an equivalent expression in their recent discussion of the 
time-independent Wigner distribution functions. 

We wish to point out that the quantum canonical transformation described here is basically different from that 
considered earlier by Kim and Wigner While the present approach deals with the transformation between 

operators {q,p) and {Q,P), their approach is about the transformation between c-numbers {q,p) and {Q,P). For the 
transformation Q — Q(q,p,t) and P = P{q,p,t), their approach yields for the kernel the expression 

<Qi,Pi,q2,P2,t) = (5[Qi - Q{q2,P2,t)]6[Pi - P{q2,P2, t)], (29) 

where Q{q,p,t) and P{q,p,t) satisfy the classical Poisson brackets relation, [Q,P]pb — [q,p]pB = 1- The kernels 
of Eq. (28) and Eq. (29) coincide with each other for the special case of a linear canonical transformation, as was 
shown by Garcia-Calderon and Moshinsky Specifically, for the case of the Wigner distribution function, they 

showed that the linear transformation for operators, Q = aq + bp and P — cq + dp, and that for c-number variables, 
Q = aq + bp and P = cq + dp, yield the same kernel k(Qi, Pi, (72,^2) = S[Qi — {aq2 + bp2)]S[Pi — {cq2 + dp2)]- 
In general cases, however, Eq. (27) and Eq. (29) give rise to different kernels. As an example, let us consider the 
unitary transformation U = e^^^^^^^ considered in Sec. II. The first-type quantum generating function has the form 



e 



iFi(quQ2)/n _ 



the kernel 



g^9{qi)/fig(^q-^ — Q2). This nonlinear canonical transformation yields for the Wigner distribution function 
/c(Qi,Pi,g2,P2) = r '^^^ I dye^l9fe-i')-9fe+y)]e¥(P2-Pi)y. (30) 

TTft j 

It is apparent that the integral of the above equation cannot generally be reduced to the (5-function of Eq. (29) except 
for some trivial cases, e.g., g =const, g = q, and g = q^. Distribution functions other than the Wigner distribution 
function do not usually allow the simple expression for the kernel in the form of Eq. (29), even if one considers a linear 
canonical transformation. 
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B. Dynamics 



In this subsection we describe how the quantum Hamihon-Jacobi theory can lead to dynamical solutions in the 
phase-space picture of quantum mechanics. For this task, we first consider the time evolution of the transformed 
distribution function in (Qi,Pi) phase space. Differentiating Eq. (25) with respect to time, we can get 



dGf 
dt 



1 



dXdYdQ2 



d 



{Q, + Y\]p\Q2-Y)t.+t{Q2 + Y\^\Q2-Y) 



d 



+t{Q2 + Y\pi^-\Q2^Y)t 
We now substitute into Eq. (31) the time evolution equations 

^^t{Q2 + Y\=t{Q2 + Y\U^, 
r) 811 - 



/(X,2r/?i)e'^('32-Qi)g-.2yPi/n_ 



where H — H{q,p,t) is the Hamiltonian governing the dynamics of the system, and obtain 
dG^ 1 



dt 



dXdYdQ2t{Q2+Y\ ( ~j-[K,p] ] \Q2~Y)tf{X,2Y/h)e 



iX{Q2~Qi)„-t2YPi/n 



(31) 

(32) 
(33) 
(34) 

(35) 



where K — K{Q, P, t) is just the transformed Hamiltonian already defined in Eq. (11). Eq. (35) should be compared 
with the following equation that governs the time evolution of the distribution function in (qi,pi) phase space, 



I 



dxdydq2{q2 + y\ --[H,p] ] \q2 -y)f{x,2y/h) 



Jx{q2-qi)p-i2ypi/h 



(36) 



We can easily see that, through the quantum canonical transformation, the role played by H is turned over to K. 

Just as the wave function has a trivial solution in the representation where the transformed Hamiltonian K(Q, P, t) 
vanishes, so does the distribution function in the corresponding phase space, as can be seen from Eq. (35). With 
the trivial solution G-^ — const., we go back to the original space via the inverse of the transformation equation 
(26) to obtain P^ {qi,pi,t). For example, for the case of the Wigner distribution fmiction the transformation can be 
accomplished by 



P^iqi,pi,t)^ J j dQ2dP2h{qi,Pl,Q2,P2,t)G^ {Q2,P2), 



where k is given in terms of the quantum principal function by 



f^^J_ I / ^y^ygi-[Si(gi+y,Q2+y,t)-S'r(9i-y,Q2-l',t)]gf [l'P2-OTi]^ 
1:1% 



(37) 



(38) 



Thus, once the quantum Hamilton- Jacobi equation is solved and the quantum principal function Si is obtained, the 
dynamics of the distribution function, as well as that of the wave function, can be determined. 



V. BOUNDARY CONDITIONS AND APPLICATIONS 

Up to this point the whole theory has been developed for the case tJ{t) = T{t)A with A taken to be arbitrary unless 
otherwise mentioned. To see how the quantum Hamilton-Jacobi theory is used to achieve the dynamical solutions of 
quantum problems, it would be sufficient, though, to consider the case of ^1 = /, the unit operator. This case was 
considered by Dirac in connection with his action principle (see Sec. 32 of Ref. [^). He showed that 5*1 defined by Eq. 
(15) equals the classical action function in the limit h 0. It should be mentioned that this particular case allows 
the quantum generating functions to attain the property that e^^^^^ is the propagator in position space and e^^^l^ the 
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propagator in momentum space. We will henceforth work on the case U{t) — T{t). The general case U{t) = T{t)A 
will be briefly treated in Appendix C. 

Before applying the theory it is necessary to provide some remarks concerning the quantum Hamilton- Jacobi 
equation (17) and its solution. First, if the Hamiltonian depends only on either q or p, we do not need to solve Eq. 
(17). Instead, since the unitary operator has the simple form U = T — e^*^*^'-'*/'' or e"'^^^-'*/'', we can obtain 5"! 
directly from Eq. (15) by calculating the matrix elements of U. For example, for a free particle, U — e~'^ t/^n^ 

is convenient to calculate e*'^^/'* = (gi|e~'^^*/^''|p2), and we get S2{qi,P2,t) = + qiP2 + i^ln2nh. Second, in 

order to solve Eq. (17), we need to impose proper boundary conditions on 5*1. Since here we are dealing with unitary 
transformations, we immediately get from the definition of Si the condition 

rfQ3e't^^(«^''33,t)-srfe,Q3,t)]M = s{qi - ga), (39) 

which follows from the calculation of the matrix elements of U{t)W{t) — I. This unitary condition ensures that 
the well-ordered operator counterpart of Si yields Hermitian operators for Q and P from Eq. (2). Mathematically, 
Eq. (17) can have several solutions, and there is an arbitrariness in the choice of the new position variable, because 
any function of the constant of integration of Eq. (17) can be a candidate for the new position variable. Not all 
the possible solutions correspond to the unitary transformations, and from the possible solutions we choose only 
those which satisfy Eq. (39) and thus give Hermitian position and momentum operators that are observables. These 
solutions correspond to the unitary transformations of the type U{t) = T{t)A. Further, from these solutions we single 
out the one that corresponds to the case A = I hy imposing the condition e^Sii<}i,Q2,t=o)/h _ _ q^-j initial 
condition. The appropriate form for iS'2 corresponding to this condition is that e^'^^^'^^'^^''^")/'' = --^l=e^'^^^^^^ . In the 
limit h 0, S2 in this equation reduces to the correct classical generating function for the identity transformation, 
S2 = qiP2- In solving the Hamilton- Jacobi equation perturbatively using Eqs. (21) and (22), in order to consistently 

satisfy the initial condition, we start with the classical Hamilton's principal function s[^^ that gives at initial time 
the relations qi — Q2 and pi — P2 from the classical c-number counterpart of Eq. (2). An arbitrary additive constant 
c to the solution of Eq. (17) that always appears in the form Si + c when we deal with a partial differential equation 
such as Eq. (19) which contains only partial derivatives of Si [Q can also be fixed by the initial condition. Depending 
whether boundary conditions can readily be expressed in a simple form, one type of the quantum generating function 
may be favored over another. The existence and uniqueness of the independent solution of Eq. (17) satisfying the 
above conditions can be guaranteed from the consideration of the equation e^Si(qi,Q2,t)/n _ (g^| J'(i)|g2), in which Si 
is just given by the matrix elements of T(t). It is clear that these matrix elements exist and are uniquely defined. 

As illustrations of the application of the theory, we consider the following two simple systems. 

Example 1. A particle under a constant force. 

As a first example, let us consider a particle moving under a constant force of magnitude a, for which the Hamiltonian 
is H = |5^/2 — aq. We start with the following classical principal function that is the solution of Eq. (21), 

^.(0) _ (gi -Q2)^ , atjqi + Q2) aH^ 

^1 - 2t + 2 ~2A- ^^^> 

Substituting 5^°^ into Eq. (22) and solving the resulting equation, we find that the first order term in h has the general 
solution 

5l^^ = ^nt + /f^-^A (41) 



2 \ t 2 



where / is an arbitrary differentiable function. To satisfy the proper boundary condition e''^i(9i,Q2,t=0)/ri _ ^^^^ 
/ and all higher order terms of Si are chosen to be zero, and the overall additive constant to be c = \ni2TTh. By 
well-ordering terms, we get the generating operator 

Si{q, Q, t) = + + + Q) „ + In ,2nht. (42) 

We can easily check that the operator form of the quantum Hamilton-Jacobi equation (13) is satisfied by the above 
generating operator. 

From Eq. (14) we obtain the wave function 

r{qi,t)^ I — ^e^[(?i~Q^)'+'^*'(?i+«^)-"'*'/i2lv-«(Q2)dQ2. (43) 
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Because ■0'^(Q2) is constant in time, we can express it in terms of the initial wave function. For the present case in 
which we use the first-type quantum generating function S'l and A = I,we have simply i/'''((?2 = Q2,t = 0) = '0'^(Q2)- 
We note that Eq. (43) is in exact agreement with the result of Feynman's path- integral approach . 

For the time evolution of the distribution function, we find from Eq. (38) the following kernel for the Wigner 
distribution function. 



ii{qi,Pi,Q2,P2,t) ^ ^^^2^ / / dYdye 



5{Q2 - 91 + Pit - af/2)5 [ P, - 2L_9l^^\ . (44) 



Substituting Eq. (44) into Eq. (37), we obtain 

F^(qi,pi,t) = ^^^(gi~Pit + atV2,Pi-ai,0), (45) 

where use has been made of the relation {qi,pi,t = 0) = {qi,pi). 

As has been mentioned, the present Hamilton-Jacobi theory also provides the solutions of the Heisenberg equations 
via the transformation relations between the two sets of canonical operators. From Eqs. (2) and (42) we can obtain 

qs = Qsit)+Ps{t)t+^t^, (46) 
Ps = Ps{t)+at. (47) 
In the Heisenberg picture, the above equations become 

qH{t)=qs+Pst+^t^, (48) 
PH{t) ^ PS + at, (49) 

which are the solutions of the Heisenberg equations. 

By setting o = 0, we can obtain the free particle solution. 
Example 2. The harmonic oscillator 

For the harmonic oscillator whose Hamiltonian is given by i? /2 + /2, the classical Hamilton-Jacobi equation 
(21) can be solved to give the classical principal function 

S'^^ ^^{qf + Ql) cot t-q^Q 2 CSC t. (50) 

With the boundary condition e''5i(«i''52,t=o)/ri ^ ^^^^ ~ Q2), Eq. (22) can be solved to give 

5^^^ = ^Insint, (51) 

(2) 

and S'l = • • • = 0. The additive constant has the form c = hii2TTh. The well-ordered generating operator is then 
written as 

- * 1 - - i 

Si{q,Q,t) = -{f + Q'^) cott - qQ csct + h-\ni2'Khsmt. (52) 

The wave function takes the form 

^"(qut)^ /^=L=e5^[(«i+^')™^*-2giQ2]^9(g2,0)dg2, (53) 
J \/i2TThsmt 

and the kernel and the distribution function are given respectively by 

K{qi,Pi,Q2,P2,t) = S{Q2 - qi cost + pi smt)S{P2 + (52 cost - qi csct), (54) 

and 

{qi,pi,t) — F^{qi cost — pi sini, qi sint + pi cost, 0). (55) 
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This equation shows that the Wigncr distribution function for the harmonic oscillator rotates clockwise in phase space. 

The quantum Hamilton- Jacobi equation for other types of generating operators can be solved by a similar technique. 
For instance, we can obtain the following solution for the second-type generating operator, 

- ' 1 - - i 

S2{q,P,t) = --{f + P^)ta.nt + qP sect + h-ln2nhcost. (56) 

The solutions of the Heisenberg equations can be obtained from Eqs. (2) and (52) (or Eqs. (3) and (56)). In the 
Heisenberg picture we have 

qnit) ^ qscost+pssint, (57) 

Ph {t) = -qs sint + ps cos t. (58) 

It should be mentioned that, even though we restricted our discussion in this section only to the case ^4 = / by 
imposing the special initial condition, it is very probable that another choice of A satisfying the quantum Hamilton- 
Jacobi equation happens to be more readily obtainable. In that case, the initial condition that is derived from 
giSi(qi,Q2fi)/h _ ^q-|^|^|g2) is of course different from that described above. As an example, for the harmonic oscillator, 
we presented a different solution for 5*1 in Appendix C where the unitary operator A corresponds to the transformation 
that interchanges the position and momentum operators. 



VI. CONCLUDING REMARKS 



We wish to give some final remarks concerning the quantum Hamilton- Jacobi theory. In this approach, the quantum 
Hamilton- Jacobi equation takes the place of the time-dependent Schrodinger equation for solving dynamical problems, 
and the quantum Hamilton's principal function that is the solution of the former equation gives the solution of 
the latter equation through Eq. (14). As mentioned in Sec. V, e^^'^l^ becomes the propagator in position space for 
the case A = I . To find the propagator, Feynman's path-integral approach divides the time difference between a 
given initial state and a final state into infinitesimal time intervals, and then lets the quantum generating function 
for the infinitesimal transformation equal the classical action function plus a proper additive constant that vanishes 
in the limit ?i ^ 0, and finally takes the sum of the infinitesimal transformations. On the other hand, the present 
approach seeks the quantum generating function that directly transforms the initial state to the final state. The 
present formalism gives also the solutions of the Heisenberg equations through the transformation relations which in 
the Heisenberg picture can be expressed as Eq. (23). In conclusion, it is clear that the present approach, which has 
its origin in Dirac's canonical transformation theory, helps better comprehend the interrelations among the existing 
different formulations of quantum mechanics. 

Finally, one more remark may be worth making as to the extent to which the quantum Hamilton-Jacobi theory 
can stretch the range of its validity. Even though our work here deals with the unitary transformation to ensure 
that the new operators become Hcrmitian, and hence observables, the main idea presented in this paper could be 
extended so as to include the non-unitary transformation that deals with non-Hermitian operators. The theory would 
then have the form of the quantum Hamilton-Jacobi equation, but it would be associated with different types of 
transformations, such as lJ(t) = T{t)B where U{t) and B are not unitary. However, it may then be necessary to pay 
particular attention and care to the completeness of the eigenstates of the new operators Q and P, for the property 
is crucial to several relations derived and has been used implicitly throughout the paper. 



APPENDIX A: PROOF OF THE EQUIVALENCE OF EQUATIONS (11) AND (12) 



Eq. (12) can be derived from Eq. (11) by considering the matrix element of the second term on the right hand side 
of Eq. (11) as follows, 

{q,\inU^\q2) = -{qi\ih^U^\q2), (Al) 
= -J dq3{qi\ih^\q,){q,\m\q2), (A2) 
= J dQs (^-ih^^{q,\U\qs)^ tiQaM, (A3) 
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where the identity UW = / is used to obtain Eq. (Al). Using the definition of the quantum generating function (1), 
we can obtain 

{q^\mj^\q2) = J dQs ^^'^l'^^^\ qimtt{Q3\q2), (A4) 
= J dQs{qi\^^^^^^\Q3)tt{Q3\q2), (A5) 

= (9i| oi k2). (A6) 

Since iTiij and g^*^'*-* have the same matrix element, we conclude that the two operators are identical. 

APPENDIX B: DERIVATION OF THE KERNEL k 

In order to find the relation between and , we first make use of the completeness of the eigenvectors of q in 
Eq. (25) and get 

Gf{QuPut) = ^J III j dXdYdQ2dq3dq^t{Q2+Y\q3){q3\p\qi){qi\Q2-Y)t 

x/(X,2F/n)e'^('5=-^i)e-'2^^i/^. (Bl) 

Changing variables with g'5 = (g3+g'4)/2 and j/ = (93—94)72, and using the relation J dy'dqeS{y—y')S{q5 — qQ)g{y',qe) = 
giUt 95) where the (5-functions are written in the forms 

Siy -y') = ^l dp3e~^'P^^y-y'y\ (B2) 

and 

%5-<Z6) = ^y^rfa;e-(«^-«=), (B3) 

we obtain 

Gf{Q,,P„t) = ^^ J J J J J J J J J dXdYdQ2dq,dydy'dp3dqedxf{X,2Y/h)t{Q2 + Y\q5+y') 

X (96 + y\p\q6 - y)(95 - y'\Q2 - y)te'^W2-0i)g-i2FPi/f.g-2ip3(j/-j/')/Re-(9«-9^). (B4) 
Next, we multiply this equation by 

where 

5{x" - x)S{y" -y) = j J dadl3e-'''^''"-''^e-^^^^y"-y^'^. (B6) 

In the resulting equation, we replace the integrations over 95 and pa, respectively, with those over 92 = 95 — a and 
P2= Pi — and then integrate over /3 and y' . We then obtain 

Gf{Qi,Put)= j j dq2dp2 J J J J j dXdYdQ2dx"dy"dat{Q2 + Y\q2 + a + y") 

j{x i^y /''') 



~2^ j j j dxdydq6{q6+y\p\qe- y)f{x,2y/h)e 
which leads immediately to Eq. (26). 



(B7) 
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APPENDIX C: SOLUTIONS OF THE QUANTUM HAMILTON-JACOBI EQUATION 



In Sec. V we considered the case A = I only for convenience, because, as can be noticed from the two examples 
of Sec. V, this case not only gives a simple relation between the initial wave fimction ijj'^{qi, 0) (distribution function 
{qi,Pi,0)) in g-representation and the constant wave function ^'^((52) (distribution function {Q2, P2)) in Q- 
representation, but also makes it easy to express the new operators {Qh,Ph) in the Heisenberg picture in terms of 
the old operators {qs,ps) in the Schrodingcr picture. In doing so, wo required the solution to satisfy the specific 
initial condition described in the first part of Sec. V. This specialization is of course not of necessity, and if this initial 
condition is discarded with the unitary condition of the transformation still retained, we have a group of general 
solutions each member of which corresponds to a specific choice of A. As an illustration, we present below another 
possible solution of the quantum Hamilton-Jacobi equation (19) for the harmonic oscillator that belongs to unitary 
transformations of the type U (t) = T{t)A, 

- - 1 - - i 

Si{q,Q,t) = --{f + Q"^) ta.nt + qQsect + n-ln2Trncost. (CI) 

Equation (CI) satisfies at initial time e''^i(9i,Q2,o)/R _ _i^^iqiQ2/h^ which equals {qi\A\q2), i.e., the matrix element 

of the unitary operator A. In this case, the operator A corresponds to the exchange transformation which generates 
the transformation relations 

P = 2^t<|p = g, (C2) 

p=-^jm^=-^. (C3) 

dQ 
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